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ABSTRACT 

We report on a detailed investigation of the high-energy 7-ray emission from 
NGC 1275, a well-known radio galaxy hosted by a giant elliptical located at the cen- 
ter of the nearby Perseus cluster. With the increased photon statistics, the center 
of the 7-ray emitting region is now measured to be separated by only 0.46arcmin 
from the nucleus of NGC 1275, well within the 95% confidence error circle with 
radius ~ 1.5arcmin. Early Fermi-L AT observations revealed a significant decade- 
timescale brightening of NGC 1275 at GeV photon energies, with a flux about seven 
times higher than the one implied by the upper limit from previous EGRET ob- 
servations. With the accumulation of one-year of Fermz'-LAT all-sky-survey expo- 
sure, we now detect flux and spectral variations of this source on month timescales, 
as reported in this paper. The average >100 MeV 7-ray spectrum of NGC 1275 
shows a possible deviation from a simple power-law shape, indicating a spectral 
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cut-off around an observed photon energy of e y = 42.2 ± 19.6 GeV, with an aver- 
age flux of F 7 = (2.31 ± 0.13) x lCT 7 phcm -2 s -1 and a power-law photon index, 
T 7 = 2.13 ± 0.02. The largest 7-ray flaring event was observed in April-May 2009 
and was accompanied by significant spectral variability above e 7 > 1 — 2 GeV. The 
7-ray activity of NGC 1275 during this flare can be described by a hysteresis behavior 
in the flux versus photon index plane. The highest energy photon associated with the 
7-ray source was detected at the very end of the observation, with the observed energy 
of e 7 = 67.4 GeV and an angular separation of about 2.4 arcmin from the nucleus. In 
this paper we present the details of the Fermi-LAT data analysis, and briefly discuss 
the implications of the observed 7-ray spectral evolution of NGC 1275 in the context 
of 7-ray blazar sources in general. 

Subject headings: galaxies: active — galaxies: jets — galaxies: individual (NGC 1275) 
— radiation mechanisms: nonthermal — gamma-rays: general 

1. Introduction 

With the successful launch of the Fermi Gamma-ray Space Telescope, we have a new op- 
portunity to study the 7-ray emission from different types of extragalactic sources — not only 
blazars, but also radio galaxies and possibly other classes of active galactic nuclei (AGN) — with 
much improved sensitivity than previously available (Abdo et al. 2010a). During the initial all-sky 
survey performed during the first four months after its launch, the Fermi Large Area Telescope 
(LAT) detected only two radio galaxies at high significance (> 10a), namely NGC 1275 (Abdo et 
al. 2009a; hereafter Paper-I) and Cen A (Abdo et al. 2009b, 2009c). More recently, the detection 
of MeV/GeV emission from yet another famous radio galaxy M 87, an established TeV source, 
was reported based on ten-months of all-sky survey Fermi-LAtl data (Abdo et al. 2009d). Yet the 
detection of NGC 1275 was particularly noteworthy because this source, unlike Cen A or M87, 
was previously undetected in 7-rays, neither by CGRO/EGRET during its ~10 years of operation, 
nor by ground-based Cherenkov Telescopes. The 7-ray flux of NGC 1275 detected by the Fermi- 
LAT was about seven times higher than the one implied by the 2a upper limit reported by EGRET, 
namely F £7> i omcV < 3.72 x 10~ 8 phcm" 2 s _1 (Reimer et al. 2003). We note that COS B data 
taken between 1975 and 1979 (Strong et al. 1982; Mayer-Has sel wander et al. 1982) showed a 
7-ray excess coincident with the position of this galaxy, although the evidence for the claimed 
high-energy source to be uniquely related to NGC 1275 is ambiguous. 

NGC 1275 is a giant elliptical galaxy located at the center of the Perseus clusteiQ. This cluster 



'The Perseus cluster (Abell 426): redshift z = 0.0179, luminosity distance cZl 
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is the brightest cluster of galaxies in the X-ray band (e.g., Bohringer et al. 1993; Fabian et al. 
2003, 2006), and as such it has been the focus of several extensive research programs over many 
years and across the entire available electromagnetic spectrum. When observed at radio wave- 
lengths, NGC 1275 hosts the exceptionally bright radio FR I radio galaxy Perseus A = 3C 84 (e.g., 
Vermeulen et al. 1994, Taylor et al. 1996, Walker et al. 2000, Asada et al. 2006). Although 
high-energy 7-rays may in general be produced within the intergalactic/interstellar medium of the 
Perseus cluster, in Paper-I, we argued that the inner radio jet of 3C 84 was the most likely source 
of the observed 7-ray photons because of the variability of the MeV/GeV flux on year/decade 
timescales implied by the EGRET and early Fermi-LAT observations. Specifically, these mea- 
surements implied the 7-ray emission region size in NGC 1275 has a radius, R < ct VSLl ~ lpc. 
Note, however, that no significant variability was detected within the four-month long Fermi-LAT 
dataset analyzed in Paper-I. Also, the LAT error circle determined from these initial data was too 
large to exclude possible contributions from other galaxies to the observed 7-ray flux (specifically, 
NGC 1273, 1274, 1277, 1278 and 1279, were still within the previously determined 95% error 
circle with R = 5.2 arcmin). 

Meanwhile, follow-up observations by the VERITAS Cherenkov telescope (Weekes et al. 
2002) put strong constraints on the very high-energy (VHE) 7-ray emission from NGC 1275 above 
100 GeV. In particular, no VHE 7-ray emission from NGC 1275 was detected by VERITAS, with 
a 99% confidence level upper limit of 2.1% of the Crab Nebula flux, corresponding to 19% of the 
power-law extrapolation of the MeV/GeV flux observed during the first four months of the Fermi- 
LAT observations (assuming the photon index, T 7 ~ 2.2; Acciari et al. 2009b). This naturally 
indicates a deviation from the pure power-law spectrum in the VHE regime, possibly compatible 
with the presence of an exponential cutoff around or below photon energies e 7 ~ 100 GeV (Ac- 
ciari et al. 2009b). The MAGIC Cherenkov telescope also recently measured upper limits for the 
VHE 7-ray emission ofNGC 1275, namely F £7> ioo GeV < (4.6-7.5) x 10~ 12 phcm~ 2 s" 1 for the 
photon indices ranging from 1.5 up to 2.5 (Aleksic et al. 2010). Thus, the implied deviation in 
the 7-ray spectrum of 3C 84 from a simple power-law form, as well as a possibility for a short- 
timescale (< month) variability of the Perseus A 7-ray flux, may now be finally addressed and 
re-examined by Fermi-LAT, due to the much improved photon statistic (especially above 10 GeV) 
after the one-year- all- sky survey. Obviously, such deep studies of NGC 1275 in the MeV/GeV pho- 
ton energy range are of the major importance for understanding the whole class of 7-ray emitting 
radio galaxies in general. 

Firmly motivated, we performed a detailed investigation of NGC 1275 in 7-rays based on the 
accumulation of one-year of Fermi-LAT all- sky survey data. In particular, we aimed to address 
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the following problems: (i) presence of short-timescale 7-ray flux variability, (ii) positional coinci- 
dence of the 7-ray emitting center with the active nucleus of NGC 1275, and (iii) spectral curvature 
and spectral evolution of NGC 1275 in the MeV/GeV photon energy range. In § 2, we describe the 
Fermi-LAT 7-ray observations and data reduction procedure. The results of the analysis are given 
in § 3, and the discussion and conclusions are presented in § 4. 



2. Fermi-LAT Observations 

The LAT instrument onboard Fermi is described in detail in Atwood et al. (2009) and refer- 
ences therein. Compared to earlier 7-ray missions, the LAT has a large effective area (~8, 000 cm 2 
on axis at 1 GeV for the event class considered here), wide energy coverage (from « 20 MeV to 
> 300 GeV), and improved angular resolution. The 68% containment angles of the reconstructed 
incoming photon direction are approximated as 6*68 ~ 0°.8 (e 7 /GeV)~ 8 below 10 GeV, giving 

~ 5°.l at 100 MeV and ~ 0°.15 at 10 GeV (update of values presented in Atwood et al. 
2009). Above 10 GeV, the improvement of angular resolution becomes relatively gentle, such that 
6> 68 ~ 0°.07 at 100 GeV. During the first year of operations, most of the telescope's time was dedi- 
cated to observing in "survey mode" where Fermi points away from the Earth and nominally rocks 
the spacecraft axis north and south from the orbital plane to enable monitoring of the entire sky on 
a time scale shorter than a day. The whole sky is surveyed every ~ 3 hours (or 2 orbits). The total 
live time included is 280.3 days (24.21 Ms). This corresponds to an absolute efficiency of 75.0 %. 
Most of the inefficiency is due to time lost during passages through the South Atlantic Anomaly 
and to readout dead time. 

The observations used here comprises all scientific data obtained between 4 August 2008 
and 13 August 2009. This time interval runs from Mission Elapsed Time (MET) 239557417 to 
271844550. We have applied the zenith angle cut of 105° to eliminate photons from the Earth's 
limb. The same zenith cut is also accounted for in the exposure calculation using the LAT science 
tooj^l GTLTCUBE. We use the "Diffuse" class events (Atwood et al. 2009), which are those recon- 
structed events having the highest probability of being photons. In the analysis presented here, we 
set the lower energy bound to a value of 100 MeV. Science Tools version v9rl5p2 and IRFs (Instru- 
mental Response Functions) P6_V3 (a model of the spatial distribution of photon events calibrated 
post-launch) were used throughout this paper. The early portion of the data here coincides with the 
early LAT observations of NGC 1275 presented in Paper-I (4 August 2008 and 5 December 2008). 



2 http://fermi. gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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2.1. Results 

Figure 1 shows the LAT 95% 7-ray localization error circles placed on the Chandra ACIS-S 
X-ray image of NGC 1275/Perseus region (ObsID 4952; exposure 164 ks). The Chandra image is 
constructed between 0.4 and 8.0 keV. The positional center of the 7-ray emission (RA= 49°. 941, 
DEC= 41°.509) is taken from the 11-month catalog (1FGLJ03 19.7+4130; Abdoetal. 2010b), and 
is only 0.46 arcmin from the NGC 1275 nucleus (RA= 49°. 951, DEC= 41°. 512), with 95% radius 
~ 1.5 arcmin (Fig. 1). Since the localization error has been substantially improved compared to 
the one quoted in Paper-I (5.2 arcmin, based on the 3-month bright source list; Abdo et al. 2009c), 
all the other field galaxies mentioned in Paper I within the 3-month LAT error (NGC 1273, 1274, 
1277, 1278, and 1279) are now outside the 11-month LAT error circle. Indeed, even the nearest 
galaxy, NGC 1274, is 2.6 arcmin apart from the positional center of the 7-ray emission with 95 
% radius of 1.5 arcmin. With the good positional association described above, we thus believe 
the 7-ray source is most likely identified as the nucleus of NGC 1275, but still needing further 
confirmation based on correlated variability with observations at other wavelengths. 

We also checked the projection of the 7-ray images as presented in Paper-I. The counts dis- 
tributions of NGC 1275 in different energy bands are consistent with the distributions for a point 
source, indicating that the diffuse extended component does not contaminate the NGC 1275/Perseus 
7-ray source, even with much improved photon statistics in the one-year dataset. Since the an- 
gular resolution of the LAT improves at high energies (Atwood et al. 2009 and see above), it 
is interesting to compare the reconstructed incoming photon directions with the position of the 
NGC 1275 nucleus around the higher energy photons that can be detected with the Fermz-LAT. 
Figure 2 shows the angular separation of e 1 > lOGeV photons from the nucleus of NGC 1275 
as a function of photon energy. The 68% and 95% angular resolutions of Fermz'-LAT are shown 
as dotted and dashed lines, respectively. The angular displacement of the highest energy photon 
detected during one-year-all-sky survey (e 7 ~ 67.4 GeV; a double circle with a number "I" in Fig- 
ure 2) from the nucleus is only 2.4 arcmin, which is well within the PSF of the LAT at this energy 
(6* 68 ~ 4.4 arcmin). Note that the chance probability for detecting Galactic and/or extragalactic 
7-ray background photons with energies at least 67.4 GeV within # 68 is less than 0.1%, hence we 
conclude that NGC 1275 is the most likely source of the discussed photon. The other three highest 
energy photons (e 7 > 30 GeV) are marked in Figure 2 as "II", "III" and "IV", and are all well 
within the 95% PSF of the Fermz'-LAT. 
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2.2. Spectral Analysis 

To study the average spectrum of NGC 1275 during the one-year observation, we use the 
standard maximum-likelihood spectral estimator provided with the LAT science tools GTLIKE. 
This fits the data to a source model, along with the models for the uniform extragalactic and 
structured Galactic backgrounds^. We use a recent Galactic diffuse model, gll_iem_v02.fit, with 
the normalization free to vary in the fit. The response function used is P6_V3_DIFFUSE. Careful 
choice of the source region is important especially for relatively faint sources. Following the 
detailed study of changing the region of interest (ROI) radius from 5° to 20° in Paper-I, we set 
r = 8° in the following analysis to minimize the contamination from the Galactic diffuse emission. 
Only one point source, corresponding to src_A in Paper-I, was found in the ROI. The coordinates 
are RA= 55°. 105, DEC= 41°. 515, with r 95 = 4.1 arcmin. We modeled the source for subsequent 
spectral analysis (also see the Fermi-LPR bright 7-ray source list; Abdo et al. 2009c). 

We first model the continuum 7-ray emission of NGC 1275 with a single power-law. The 
extragalactic background is assumed to have a power-law spectrum as well, with the spectral index 
and the normalization free to vary in the fit. From an unbinned GTLIKE analysis, the best fit 
power-law parameters for NGC 1275 are: 



^ = (2.61 ±0.14) x 10- 9 ( \\ - )-2-i3±o.o2 phcm-^-iMeV" 1 , (1) 
cfe 7 100 MeV 



or 



i^iooMeV = (2.31 ±0.13) x 10~ 7 phcm-V 1 , (2) 

with only statistical errors taken into account. Systematic errors for the LAT are still under investi- 
gation , but the estimated systematic uncertainty on the flux is 10 % at 100 MeV, 5 % at 500 MeV, 
and 20 % at 10 GeV, respectively (Abdo et al. 2010b). The results are consistent with those pre- 
sented in Paper-I, even though the uncertainties are now much smaller due to the improved photon 
statistics. 

The predicted photon counts from NGC 1275 in the ROI are V prc d = 3187.0 and the test 
statistic (defined as TS = 2 [logL — logL ], where L and L are the likelihood when the source 
is included or not, respectively) is TS = 4039.8 above £ 7 = 100 MeV, corresponding to a 64cx 
detection. For the Galactic diffuse background, the normalization is 1.11 ± 0.01 and N pre d = 
44835.3. The power-law photon index of the extragalactic background is T 7 = 2.31 ± 0.03 with 



http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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Np Te d = 9081.7. Figure 3 shows the LAT spectrum of NGC 1275 obtained by separately running 
GTLIKE for 10 energy bands; 100 - 200 MeV, 200 - 400 MeV, 400 - 800 MeV, 800 MeV -1.6 GeV, 
1.6 - 3.2 GeV, 3.2 - 6.4 GeV, 6.4 - 12.8 GeV, 12.8 - 25.6 GeV, 25.6 - 51.2 GeV, and 51.2 - 
102.4 GeV, where the dashed line shows the best fit power-law function for the NGC 1275 data 
given in equation (1). For the highest energy bin (51.2 — 102.4 GeV), we plot a 2 o upper limit 
since this bin includes only one photon (see Figure. 2) and is ignored in the subsequent statistical 
analysis. 

Note an indication of a deviation of the model with respect to the data above e 7 = 20 GeV. 
Indeed, a x 2 fit of the power-law model to the datc@ gives a relatively poor fit with x 2 = 14.0 for 7 
degrees of freedom (d.o.f.), where its probability is P(x 2 ) — 0.05. Instead, an alternative fit with 
a cutoff power-law function with the form dN/de^ oc e 7 7 x exp [— e^/e c ] gives a much better 
representation of the data with an improved x 2 = 6.7 for 7d.o.f., corresponding to P(x 2 ) = 0.34. 
We therefore retried the GTLIKE fit assuming a power-law with an exponential cutoff function, 
and obtained e c = 42.2 ± 19.6 GeV with T 7 = 2.07 ± 0.03. Moreover, we apply a likelihood 
ratio test between a simple power-law and a cutoff power-law function. The test statistic (D) is 
twice the difference in these log-likelihoods, which gives D = 2.9 for our case. Note that the 
probability distribution of the test statistic can be approximated by a x 2 distribution with 1 d.o.f., 
corresponding to different d.o.f. between two functions. We obtain P(x 2 ) < 0.08, which again 
indicates a deviation from a simple power-law function although this is currently inconclusive. The 
best-fit cutoff power-law function is shown as dotted line in Figure 3. 



2.3. Temporal Variability 

Next, we investigated the 7-ray flux variations of NGC 1275 from 4 August 2008 to 13 August 
2009. To this end, we constructed spectra with a time resolution of 14 days and fit each spectrum 
with a power-law model just for simplicity. The ROI radius (r = 8°), the energy range (e 7 > 
100 MeV), and other screening conditions are the same as described above. Since no variability 
is expected for the underlying background diffuse emission, we fix the best-fit parameters to the 
average values determined from the one-year integrated spectrum for the Galactic and extragalactic 
background components. We first fix the spectral index to the best fit value over the full interval, 
T 7 = 2.13, to minimize uncertainties in the flux estimates. In this case, variability is highly 
significant with x 2 — 157.7 for 26 d.o.f., where P(x 2 ) < 10~ 6 . Indeed, the 7-ray flux varies 



Note that the maximum likelihood itself does not provide any information about the quality of a fit for an assumed 
model. We therefore perform x 2 fitting to the resultant 10 band LAT spectrum as described above, to give a convenient 
estimate of the goodness of the fit. 
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even within a few month timescale. Due to the limited photon statistics, however, it is difficult 
to investigate any shorter timescale variability of NGC 1275. We also investigate the spectral 
evolution of NGC 1275, with the 7-ray photon index free to vary. Figure 4 shows variations of the 
flux (e 7 > 100 MeV: upper panel) and photon index {lower panel) versus time. A flaring event is 
seen around T = 252 — 294 days after 2008 August 4th (epoch "B" in Figure 4), corresponding 
to the epoch April-May 2009. A doubling timescale for the flaring period cannot be accurately 
measured, but ~20 days. A x 2 fit to a constant gives 35.6 and 43.0 for 26d.o.f., corresponding to 
the probability P(x 2 ) = 0.10 and 0.02, respectively, for the flux and the photon index variations. 
Interestingly, the spectral variation is predominantly due to the difference between the pre- and 
post-flare period, denoted in Figure 4 as epochs "A" and "C". In fact, if we test the variability 
separately for the pre- and post-flare periods (A and C), the significance of the variability decreases 
substantially, such that a constant fit provides P{x 2 ) > 0.21 or 0.29 for both the flux and the photon 
index changes. During the flare (epoch B), the photon index changes from T 7 ~ 2.2 to T 7 < 2.0. 
Note that, after the flare (epoch C), the flux drops to its original level but the spectrum remains 
relatively flat, with T 7 ~ 2, persisting for more than a few months. It should also be noted that 
three out of the four highest energy photons (denoted as I, II, and III in Figure 2) were indeed 
detected in the post-flare period, when the spectrum appears the flattest. 



3. Discussion and Conclusion 

In the previous sections, we reported on the analysis of the 7-ray emission from NGC 1275 
observed with Fermi-LAI during its one-year- all- sky survey. We showed that with the increased 
photon statistics, the positional center of the 7-ray emission is now much closer to the NGC 1275 
nucleus as compared to that reported in Paper-I. In addition, we have shown that the average 7- 
ray spectrum of NGC 1275 reveals a significant deviation from a simple power-law above photon 
energies e y ~ 1 — 2 GeV. That is, the observed Fermi-LAI spectrum is best fitted by a power-law 
function (r 7 ~ 2.1) with an exponential cutoff at the break photon energy e c = 42.2 ± 19.6 GeV. 
Finally, we argued that significant flux and spectral changes of NGC 1275 are detected with Fermi- 
LAI on a timescale of a few months, although the possibility for even shorter variability remains 
uncertain. 

We also reported the detection of an interesting spectral evolution, consisting of a persistent 
(over more than a few months) spectral hardening (from T 7 ~ 2.2 to T 7 ~ 2) after the largest 
flaring event observed in April-May 2009. During this flat-spectrum/low-flux-level epoch the 
highest energy photon (s 1 ~ 67.4 GeV) was detected from the direction of NGC 1275. All these 
new findings basically support the idea put forward in Paper-I that the observed 7-ray emission 
from the Perseus system originates in the (sub) pc-scale radio jet of NGC 1275, and is therefore 
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most likely analogous to high-energy emission observed in blazars. In fact, as shown in Paper-I, 
the overall vF v spectral energy distribution (SED) of NGC 1275 constructed with multi-frequency 
radio to 7-ray data shows a close similarity to the "two-bump" SEDs of so-called low-frequency 
peaked BL Lac objects (hereafter LBLs). Till now, only a few LBLs have been detected at TeV 
photon energies - BL Lac (Albert et al. 2007), 3C66A (Acciari et al. 2009a), S5 0716+714 
(Anderhub et al. 2009), and also W Comae from an IBL class (Acciari et al. 2008; 2009d) - and 
more similar discoveries are expected. Hence, NGC 1275 itself has been suggested to be a potential 
TeV source as well, thus motivating deep VERITAS (Acciari et al. 2009b) and MAGIC (Aleksic 
et al. 2010) observations. These observations so far resulted only in upper limits to the VHE 7-ray 
emission of the studied region. One should note, however, that both the low- and high-energy 
peaks of NGC 1275 (in the vF v representations) are located at substantially lower frequencies than 
those of typical LBLs. Indeed, the low-energy (synchrotron) emission components of 3C 66A and 
BL Lac peak around 10 13_15 Hz, while it is around 10 12 Hz in the case of NGC 1275 (Paper-I). 
Correspondingly, the high-energy emission component of 3C 66A (but not necessarily of BL Lac) 
peaks at higher photon energies than that observed in NGC 1275. Within such an interpretation, 
this is consistent with the softer MeV/GeV spectrum in NGC 1275 (1% = 2.13 ± 0.02) compared 
to 3C66A (r 7 = 1.97 ± 0.04; Abdo et al. 2009b) as measured by the Fermi-LAT. Note also 
that the X-ray spectrum of NGC 1275 is clearly "rising" in the vF v representation (r x = 1.6—1.7; 
Balmaverde et al. 2006, Ajello et al. 2009), indicating that it is dominated by the low-energy 
portion of the IC emission, whilst the X-ray spectra of LBLs are typically very flat, suggesting a 
transition between the synchrotron and IC components (I\ — 2; e.g., Ghisellini et al. 1998). Such 
observational differences may indicate that, unlike in the case of LBLs, the high-energy spectrum 
of NGC 1275 does not extend up to TeV photon energies. 

An interesting comparison can be made to another nearby radio galaxy detected by Fermi- 
LAT, namely M 87 which resembles NGC 1275 in many respects, and is an established TeV source 
(Abdo et al. 2009d). While M87 is closer to us than NGC 1275 (d L = 16Mpc versus d L = 
75.3 Mpc), the GeV flux of M87 is much lower than that of NGC 1275 (F £7>100 McV = (2.45 ± 
0.63) x 10- 8 phcm- 2 s _1 versus (2.31 ± 0.13) x 10~ 7 phcm- 2 s" 1 ), with similar LAT measured 
spectra (r 7 = 2.26 ± 0.13 versus 2.13 ± 0.02). Note also that both radio galaxies are located at 
the centers of rich clusters, that the synchrotron emission components in both sources peak in the 
far infrared (~ 10 12 — 10 13 Hz), and the estimated jet powers are similar, L } ~ 10 44 erg s _1 (Owen 
et al. 2002, Dunn & Fabian 2004). Thus, it may be surprising that only one of these has so far 
been detected at TeV photon energies. The detailed analysis of the spectral evolution of NGC 1275 
within the Fermi-LAT range reported in this paper may provide a viable explanation for such a 
behavior. In particular, as already emphasized above, we found that the epochs characterized 
by the flattest GeV continuum of this source, as well as the arrival times of the highest energy 
photons from the direction of NGC 1275, do not coincide with the epochs of the highest photon 
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flux above 100 Me V. This is clearly illustrated in Figure 5, which shows a correlation between 
F £7>100 McV an d the photon indices emerging from the power-law fits. Because no significant flux 
or spectral changes were observed during the pre- and post-flare epochs A (0 — 252 day) and C 
(294 — 374 day) (see Figure 4), the average fluxes and photon indices for these time periods are 
reported. Note that the average fluxes of these pre- and post-flare epochs are comparable, whilst the 
corresponding photon indices differ significantly by Ar 7 ~ 0.2. Moreover, the observed spectral 
evolution during one year of the Fermi-LAT exposure reveals a hysteresis-like character, more 
clearly seen for the flaring period (epoch B, with a time bin of 2 weeks), followed by a gradual 
flattening in the subsequent decay phase. Further insights into the spectral evolution of NGC 1275 
within the Fermi-LAT photon energy range are provided by Figure 6, which shows the two SEDs 
for pre- and post-flare epochs A and C. Here the dotted line corresponds to the best "power-law 
with an exponential cutoff" fit function determined from an average 7-ray spectrum, as given 
in Figure 3. Interestingly, the difference between the two SEDs consists of an excess at photon 
energies £ 7 ~ 1 — 2 GeV, with the low-energy 7-ray flux remaining essentially unchanged between 
the two epochs. This implies that (1) the 7-ray variability in NGC 1275 (and possibly other radio 
galaxies) may be restricted to > GeV photon energies, and that (2) the position of the peak in the 
high-energy spectral component (in the vF v representation) may change substantially even within 
the same object with no accompanying significant flux changes. Note in this context that both the 
VERITAS and MAGIC non-detections were obtained during the pre-flare epoch A (164 — 206 day; 
see Acciari et al. 2009b; Aleksic et al. 2010). Hence, the emerging conclusion is that it is not 
the total flux above 100 MeV which should play a major role in triggering TeV observations of 
steep-spectrum Fermi-LAT sources, but instead it is the flux and photon index determined at higher 
photon energies (> GeV). 

On the theoretical side, this conclusion could be possibly justified by noting that after a new 
episode of injection of freshly accelerated electrons into the emission zone (e.g., a downstream 
region of a shock), higher energy electrons may lag behind the lower-energy electrons. In such a 
case, as discussed previously in the context of blazar modeling (e.g., Kirk et al. 1998; Sato et al. 
2008), an "anti-clockwise" hysteresis in the flux versus photon index plane may arise, similar to 
what we observed in NGC 1275 at 7-ray photon energies. 

Here, we comment on the question if NGC 1275 - being a representative example of a low- 
power radio galaxy - may be considered as a misaligned blazar, most likely of the LBL type. 
Assuming a homogeneous jet model, one should expect its jet Doppler factor 5 = rr 1 (1 — 

/3j cos = iy 1 (1 — y^l — iy 2 cos 6)~ l ~ 1 — 2, for the typically expected jet viewing angle 
of 6 ~ 20° — 30° and jet bulk Lorentz factors Tj ~ 10. Indeed, modeling of the broad-band 
emission of LBLs (beamed counterparts of radio galaxies such as NGC 1275 by assumption) re- 
quires 5 ~ Tj ~ 10. Thus, if the only difference between radio galaxies and blazars is due to the 
viewing angle, this should manifest in (i) different observed positions of the spectral peaks in the 
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vF v representation {y oc 5), (ii) different observed variability patterns (t var oc assuming the 
emission region is a moving source), and finally in (iii) different observed luminosities (L obs oc 5 4 
for a moving blob case, or oc 5 3 /Tj for the steady jet (see Sikora et al. 1997). 

Indeed, both the low- and high-energy vF v peaks of NGC 1275 are located at substantially 
lower frequencies than those of typical LBLs. Taking the difference of beaming factors into ac- 
count (i.e., 5e =0 /5e=2o° ~ 10), the broad-band SED of the NGC 1275 would be similar to the 
SEDs of blazars such as BL Lacertae or 3C 66A. Note however that the position of the high-energy 
spectral peak may change substantially in a single object even for comparable flux levels, at least 
in NGC 1275, so the diagnostics related to the location of the spectral peaks may not be very 
conclusive. 

Variability as short as day timescales is often observed in LBLs. For example, during the 
historical flare of BL Lacertae in 1997, correlated 7-ray and optical flares were observed, with the 
7-ray flux increasing by a factor of 2.5 within a day (Bloom et al. 1997). Similarly, daily variability 
has been discovered in both Fermi-LAI and VERITAS observations of 3C 66A (Reyes et al. 2009). 
Hence, in the simple unification scheme outlined above, we should expect NGC 1275 to vary in 
7-rays on timescales of several days. This, however, is difficult to test even with the excellent 
sensitivity of the Fermi-LAI instrument, due to the limited photon statistics for weekly time bins. 
The analysis presented in this paper gives instead only a robust upper limit t var < a few months. 
However, note that day-timescale variability has been detected at TeV photon energy range for 
the M87 radio galaxy (Acciari et al. 2009c). Thus, more frequent monitoring of NGC 1275 by 
ground-based Cherenkov Telescopes would be valuable in this respect. 

The observed photon flux of -F £7 >iooMeV = (2.19 ± 0.13) x 10~ 7 phcm~ 2 s _1 implies an ob- 
served (isotropic) 7-ray luminosity of NGC 1275, L 7 ~ Andl (r 7 — 1) F £i>£0 (e / e o y~ ri de ~ 
10 44 ergs _1 , for T 7 ~ 2.0 — 2.2 and e c ~ 42GeV. This is already comparable to the typical 
observed 7-ray luminosities of LBLs, which range between L 7 ~ 10 44 — 10 46 erg s -1 . Note in par- 
ticular that in the framework of the simple unification scheme discussed above, the beamed analog 
of NGC 1275 would then be characterized by the observed 7-ray luminosity larger by a factor of 

[5 e=0 /5o =2 o°] 4 ~ 3 x 10 4 (or at least [5 e=Q /5 e =2o°f ~ 3 x 1q3 ) tnan tnis > i- e - L ~i > 10 47 erg s" 1 . 
Such luminosities are not expected for LBL-type blazars (see Abdo et al. 2010a). On the other 
hand, the observed 7-ray luminosity of NGC 1275 is not energetically problematic, since the total 
emitted 7-ray power in this source seems rather moderate as long as the emitting plasma moves 
with highly relativistic bulk velocities, namely L 7 e m — (Oj/47r) L 7 ~ L 7 /4r| < 10 42 ergs _1 
for Tj ~ 10, where Qj ~ ndf is the solid angle defined by the jet opening angle 0- p for which 
we assumed 0j ~ Such a relatively small emitted power would constitute less than 1% of 

the total kinetic power of the NGC 1275 jet, estimated by Dunn & Fabian (2004) to be roughly 
Lj ~ (0.3 — 1.3) x 10 44 ergs _1 . Yet the problem of an unexpectedly large observed 7-ray lumi- 
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nosity of the beamed analog of NGC 1275 remains, and poses a serious challenge to the simplest 
version of the AGN unification scheme. In this context, a viable explanation for this problem 
would be to postulate that the high-energy emission observed from "misaligned" blazars such as 
NGC 1275 (or M 87) is dominated not by a jet "spine" characterized by large bulk Lorentz factors 
(Tj ~ 10 as is the case in bona-fide blazars), but by the slower jet boundary layers (Tj ~ few) as 
discussed by several authors (e.g., Celotti et al. 2001, Stawarz & Ostrowski 2002, Ghisellini et 
al. 2005). Alternatively, one may propose that the 7-ray emission observed from radio galaxies 
is not produced within the "proper" blazar emission zone, but at larger distances from the active 
center characterized by slower bulk velocities (say, Tj ~ few) of the emitting plasma. Yet an- 
other possibility may be that the inner jets in NGC 1275 (and also in similar objects) are in general 
intrinsically less relativistic than the ones in bona-fide blazars; this would be consistent with the 
conclusions of Lister & Marscher (1997), who argue that radio-loud AGN with nuclear jets char- 
acterized by Tj > 10 must be rather rare among the general population. Whichever scenario is 
correct, the Fermi results seem to indicate that low-power radio galaxies are most likely not simple 
off-axis analogs of BL Lac objects in terms of their 7-ray properties. 
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Fig. 1. — The 95 % LAT 7-ray localization error circles placed on the Chandra ACIS-S image 
constructed between 0.4 and 8.0 keV. The positional center of the 7-ray emission, marked as the 
white "X", is only 0.46 arcmin from the position of the NGC 1275 nucleus, with 95% radii of r 95 
= 5.2 arcmin for Paper-I and 1.5 arcmin (more accurately, 1.56 x 1.38 arcmin; Abdo et al. 2010b) 
for 11 months, respectively. 
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Fig. 2. — Angular separation of high energy (e 7 > 10 GeV) photons from the nucleus of NGC 1275 
as a function of photon energy. The 68% and 95% angular resolutions of Fermi-LAI are shown 
as dotted and dashed lines, respectively. These PSF profiles are derived for P6_V3_Diffuse and 
have been averaged over the acceptance of the LAT. The highest energy photon detected during the 
one-year- all- sky survey is £ 7 = 67.4 GeV (double circle denoted as "I"), whose angular separation 
from NGC 1275 is only 2.4 arcmin. Arrival times of the four highest energy photons ("I-IV"), are 
indicated in Figure 4. 
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Fig. 3. — The average LAT >100 MeV spectrum of NGC 1275 derived from one year accumula- 
tion of data (filled circles; this work), as compared with that determined from the initial 4-month 
dataset (open squares; Paper I). The dashed line shows the power-law function for the one-year 
data determined from the GTLIKE. The dotted line represents the best- fit cutoff power-law with 
e c = 42.2 GeV as described in the text. 
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Fig. 4. — Temporal variation of 7-ray flux and spectral index over the period August 2008 - August 
2009. The time (in days) is measured from the start of the Fermi observation, i.e., 2008 August 
4, 15:43:37 UT. Upper panel: changes in the £ 7 > 100 MeV flux. Lower panel: changes in the 
power-law photon index. We have divided the analyzed time window into epochs A (before the 
flare), B (during the flare), and C (after the flare). The arrival times of the highest energy photons 
I-IV (Figure 2) are indicated as arrows. Background diffuse emission (both Galactic and extra- 
galactic) is fixed at the best-fit parameters determined from an average spectral fit as given in the 
text, and only statistical errors are shown. 
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Fig. 5. — Spectral evolution of NGC 1275 in the flux-photon index representation. Labels "A" and 
"C" denote the average fluxes and photon indices determined for the period before (0—252 day) and 
after the flare (294 — 374 day), respectively (see Figure 4). Label "B" denotes the spectral evolution 
during the flare with a time bin of 2 weeks. Note that the photon index changes significantly before 
and after the flare, while the flux levels are comparable. 
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Fig. 6. — The 7-ray spectral energy distribution of NGC 1275 before (crosses) and after the flare 
(open circles), corresponding to epochs A and C in Figure 4, respectively. Note that the excess 
high-energy 7-ray emission only appears above £ 7 ~ 1 — 2 GeV. The flux upper limits in the TeV 
range resulting from the VERITAS and MAGIC observations correspond to the pre-flare epoch A. 
The dotted line represents the best- fit cutoff power-law fit as given in Figure 3. 



